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ABSTRACT 


This  NISTIR  serves  as  the  user’s  manual  for  Version  1.0  of  the  Virtual  Cement  and  Concrete  Test- 
ing Laboratory  (VCCTL).  Using  the  VCCTL,  a user  may  create  starting  microstruct.ures  of  cement 
(gypsum,  fly  ash,  etc.)  particles  in  water,  hydrate  the  microstructures  under  controlled  curing  con- 
ditions, leach  the  portlandite  and  remaining  cement  clinker  phases  from  hydrated  microstructures, 
view  images  and  animations  of  3-D  microstructures,  and/or  evaluate  a variety  of  properties  of  hy- 
drated/leached microstructures  for  direct  comparison  to  experimental  data.  The  VCCTL  includes 
a web-based  interface  that  provides  a menu-driven  system  of  forms  for  obtaining  user  input  and 
returning  results.  Two  example  applications  of  using  the  VCCTL  to  study  specific  scenarios  of 
interest  to  the  cement-based  materials  research  community  are  also  presented. 

Keywords:  Building  technology,  cement  hydration,  computer  modeling,  concrete  testing,  mi- 
crostructure, simulation,  virtual  laboratory. 
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1 Introduction 


This  user’s  manual  provides  guidance  for  using  Version  1.0  of  the  NIST  Virtual  Cement  and  Con- 
crete Testing  Laboratory  (VCCTL).  This  virtual  laboratory  consists  of  a web-based  interface  and 
underlying  computer  models  and  programs  that  allow  the  user  to  create,  hydrate,  degrade,  and  eval- 
uate 3-D  cement-based  microstructures  from  their  desktop  computer.  Thus,  the  VCCTL  should 
serve  the  cement  and  concrete  research  and  development  community,  ultimately  reducing  the  num- 
ber of  costly  physical  concrete  tests  needed,  while  also  allowing  the  rapid  analysis  of  a wide  variety 
of  “what  if”  type  scenarios. 

The  main  underlying  program  for  the  VCCTL  (CEMHYD3D)  has  been  described  previously  [1], 
as  well  as  numerous  applications  based  on  using  the  underlying  models  [2-7].  This  user’s  manual 
therefore  provides  a detailed  description  of  the  web-based  user  interface  and  a few  examples  of 
using  the  VCCTL  to  investigate  specific  problems  of  interest  to  the  cement  and  concrete  research 
community.  The  VCCTL  is  available  online  at  http://vcctl.cbt.nist.gov.  It  is  suggested  that  the 
user  be  connected  to  this  web  site  as  they  read  this  user’s  manual  in  order  to  help  clarify  the  forms, 
the  inputs  needed,  and  the  results  produced  for  the  various  menu  selections. 

DISCLAIMER 

This  software  was  developed  at  the  National  Institute  of  Standards  and  Technology  by  employees 
of  the  Federal  Government  in  the  course  of  their  official  duties.  Pursuant  to  title  17  Section  105 
of  the  United  States  Code  this  software  is  not  subject  to  copyright  protection  and  is  in  the  public 
domain.  The  VCCTL  is  an  experimental  system.  NIST  assumes  no  responsibility  whatsoever  for 
its  use  by  other  parties,  and  makes  no  guarantees,  expressed  or  implied,  about  its  quality,  reliability, 
or  any  other  characteristic.  We  would  appreciate  acknowledgement  if  the  software  is  used. 

The  U.S.  Department  of  Commerce  makes  no  warranty,  expressed  or  implied,  to  users  of  the 
VCCTL  and  associated  computer  programs,  and  accepts  no  responsibility  for  its  use.  Users  of  the 
VCCTL  assume  sole  responsibility  under  Federal  law  for  determining  the  appropriateness  of  its  use 
in  any  particular  application;  for  any  conclusions  drawn  from  the  results  of  its  use;  and  for  any 
actions  taken  or  not  taken  as  a result  of  analyses  performed  using  these  tools. 

Users  are  warned  that  the  VCCTL  is  intended  for  use  only  by  those  competent  in  the  field  of 
cement-based  materials  and  is  intended  only  to  supplement  the  informed  judgment  of  the  qualified 
user.  The  software  package  contains  computer  models  which  may  or  may  not  have  predictive  value 
when  applied  to  a specific  set  of  factual  circumstances.  Lack  of  accurate  predictions  by  the  models 
could  lead  to  erroneous  conclusions  with  regard  to  materials  selection  and  design.  All  results  should 
be  evaluated  by  an  informed  user. 


INTENT  AND  USE 

The  algorithms,  procedures,  and  computer  programs  described  in  this  report  constitute  a pro- 
totype system  for  a virtual  laboratory  for  the  testing  of  cement  and  concrete.  They  have  been 
compiled  from  the  best  knowledge  and  understanding  currently  available,  but  have  important  lim- 
itations that  must  be  understood  and  considered  by  the  user.  The  VCCTL  system  is  intended  for 
use  by  persons  competent  in  the  field  of  cement-based  materials  and  with  some  familiarity  with 
computers.  It  is  intended  as  an  aid  in  the  materials  selection,  optimization,  and  design  process. 

2 Descriptions  of  Menu  Selections 

Upon  entering  the  VCCTL  system,  the  main  menu  shown  in  Figure  1 is  displayed.  Each  of  the 
possible  main  menu  selections  will  be  discussed  in  the  subsections  which  follow. 
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Main  Menu  for  Virtual  Cement  and  Concrete 

Testing  Laboratory 


Please  select  from  the  following  menu: 

0)  Access  the  User’s  Guide  for  the  VCCTL 

1)  Enter  the  cement  images  database  system 

2)  Create  a particle  size  distribution  file 

3)  Generate  an  initial  3-D  microstructure 

4)  View  a 2-D  slice  from  any  3-D  microstructure 

5)  Generate  an  animation  from  any  3-D  microstructure 

6)  View  an  animation  from  any  3-D  microstructure 

7)  Distribute  all  phases  amongst  cement  particles 

8)  Filter  a 3-D  microstructure  to  distribute  only  two  phases  amongst  cement 
particles 

9)  Compute  spatial  statistics  for  a 3-D  microstructure 

10)  Sinter  a 3-D  microstructure  to  adjust  the  hydraulic  radius  of  one  phase 

11)  Distribute  fly  ash  phases  amongst  fly  ash  particles 

12)  Input  a temperature  history  (curing)  file 

13)  Hydrate  a 3-D  microstructure 

14)  Measure  phase  fractions  vs.  distance  from  aggregate 

15)  Leach  CH/cement  from  a hydrated  microstructure 

16)  Evaluate  percolation  properties  of  a hydrated/leached  microstructure 

17)  Plot  predicted  properties  (vs.  experimental  data) 

18)  Predict  chloride  ion  diffusivity  of  a concrete  based  on  mixture  parameters 

Figure  1:  Main  menu  for  the  Virtual  Cement  and  Concrete  Testing  Laboratory. 
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2.1  Menu  Selection  0)  Access  the  User’s  Guide  for  the  VCCTL 


Menu  Selection  0 simply  provides  a link  to  the  HTML  version  of  this  user’s  guide  available  online 
as  part  of  the  VCCTL  system. 

2.2  Menu  Selection  1)  Enter  the  cement  images  database  system 

For  users  of  the  VCCTL  who  do  not  have  access  to  a scanning  electron  microscope  (SEM)  with 
X-ray  imaging  capabilities,  a database  of  cement  images  has  been  created.  Currently,  the  database 
contains  2-D  images,  quantitative  image  analysis  residts,  and  measured  particle  size  distributions 
(PSDs)  for  17  different  cements  produced  in  four  different  countries,  as  shown  in  Figure  2.  The 
characterization  process  has  been  described  in  detail  for  two  of  the  Cement  and  Concrete  Reference 
Laboratory  cements  in  a previously  published  report  [6].  As  new  cements  are  processed  at  NIST, 
they  will  be  added  to  the  existing  database.  Typically,  the  PSD  and  quantitiative  phase  analysis 
information  will  be  used  as  input  to  construct  a starting  3-D  microstructure  for  a specific  cement. 

Upon  entering  the  cement  images  database,  the  user  will  first  be  allowed  to  select  a specific 
cement  of  interest  using  a fill-in  form.  Once  this  form  is  submitted,  a 2-D  color  processed  SEM/X- 
ray  image  will  be  returned  along  with  the  measured  PSD  for  the  chosen  cement,  as  shown  in  Figures 
3 and  4 for  Cement  135,  issued  by  the  Cement  and  Concrete  Reference  Laboratory  in  January  of 
2000.  The  PSD  is  reported  on  a particle  diameter  basis,  with  both  mass  and  number  fractions 
reported.  The  number  fractions  have  been  obtained  from  the  mass  fractions  based  on  the  number 
of  pixel  elements  present  in  digitized  3-D  particles  of  each  specific  diameter  [1].  For  small  diameters 
(3  pixels  and  5 pixels),  it  should  be  noted  that  the  volume  of  the  digitized  sphere  is  significantly 
greater  than  an  equivalent  diameter  continuum  sphere. 

Each  cement  page  also  contains  two  links  to  further  information.  The  About  this  cement 
link  directs  the  user  to  a secondary  page  that  contains  the  quantitative  image  analysis  results 
for  this  cement,  including  the  names  of  the  computed  correlation  files  which  can  be  used  to  cre- 
ate “equivalent”  3-D  microstructures.  Normally,  this  page  will  be  printed  by  the  user  for  later 
use  with  Menu  Selection  7 to  distribute  the  cement  phases  amongst  the  3-D  particles.  The 
Access  the  ftp  site  for  these  files  link  directs  the  user  to  an  anonymous  ftp  site  from  which  the 
correlation  and  PSD  files  can  be  downloaded  for  personal  use.  Optionally,  the  measured  PSD  pro- 
vided on  the  cement  page  can  be  printed  and  typed  in  by  the  user  as  input  into  Menu  Selection  2) 
Create  a particle  size  distribution  file. 
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Cement  Images  and  PSDs  Available  for  Microstructure  Modeling 


Please  select  one  of  the  portland  cements  below  to  view  image  and  PSD  info 
(then  press  the  submit  button  at  the  bottom  of  the  page). 

CCRL  Cement  1 15  (Jan.  1995) 

CCRL  Cement  116  (Jan.  1995) 

CCRL  Cement  133  (June  1999) 

CCRL  Cement  135  (Jan.  2000) 

CCRL  Cement  136  (Jan.  2000) 


Montalieu  Cement  (French) 

South  African  Cement 
Dyckerhoff  Cement  654  rrr/kg  (Germany) 
Dyckerhoff  Cement  370  rrr/kg  (Germany) 
Dyckerhoff  Cement  197  irr/kg  (Germany) 

'y 

Dyckerhoff  Cement  608  rrr/kg  (Germany) 
Dyckerhoff  Cement  357  rrr/kg  (Germany) 
Dyckerhoff  Cement  248  irr/kg  (Germany) 
Submit  form  for  response 


>y 

Dyckerhoff  Cement  541  nr/kg  (Germany) 

'y 

Dyckerhoff  Cement  279  rrr/kg  (Germany) 

'y 

Dyckerhoff  Cement  468  rrr/kg  (Germany) 

'y 

Dyckerhoff  Cement  281  rrr/kg  (Germany) 


Return  to  Homepage  for  Partnership  for  High  Performance  Concrete  Technology 


Figure  2:  Form  to  select  a cement  from  the  cement  images  database. 
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cementl35 


Figure  3:  First  half  of  page  from  cements  image  database  for  Cement  135. 
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Discretized  Cement  Particle  Size  Distribution 


Diameter_(um) 

Wt._frac. 

Part_vol 

Number_frac 

1.0 

0.177775 

1 

0.956205 

3.0 

0.121252 

19 

0.034325 

5.0 

0.087078 

81 

0.005782 

7.0 

0.067035 

179 

0.002014 

9.0 

0.055533 

389 

0.000768 

11.0 

0.047494 

739 

0.000346 

13.0 

0.042191 

1189 

0.000191 

15.0 

0.038665 

1791 

0.000116 

17.0 

0.036427 

2553 

0.000077 

19.0 

0.032629 

3695 

0.000047 

21.0 

0.031640 

4945 

0.000034 

23.0 

0.028675 

6403 

0.000024 

25.0 

0.027466 

8217 

0.000018 

27.0 

0.023840 

10395 

0.000012 

29.0 

0.023840 

12893 

0.000010 

31.0 

0.028647 

15515 

0.000010 

35.0 

0.041058 

22575 

0.000010 

41.0 

0.035907 

36137 

0.000005 

47.0 

0.031764 

54435 

0.000003 

61.0 

0.016722 

119009 

0.000001 

73.0 

0.004360 

203965 

0.000000 

Access  ftp  site  for  these  files 


Figure  4:  Second  half  of  page  from  cements  image  database  for  Cement  135. 
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2.3  Menu  Selection  2)  Create  a particle  size  distribution  file 

This  menu  selection  is  used  to  create  and  store  a PSD  file  which  specifies  the  number  of  cement 
particles  of  each  diameter  to  be  placed  in  a starting  3-D  microstructure.  The  output  file  created  by 
this  selection  can  be  specified  as  input  to  Menu  Selection  3)  Generate  an  initial  3-D  microstructure. 

The  form  (shown  in  Figure  5)  to  create  a PSD  file  consists  of  a table  listing  the  diameter,  mass 
fraction,  and  number  of  particles,  followed  by  two  entries  to  specify  the  total  number  of  particles 
or  the  water- to-cement  ratio  (w/c),  respectively.  The  user  may  edit  the  mass  fractions  to  match 
those  obtained  from  one  of  the  cements  in  the  cements  image  database  as  discussed  above.  Then, 
they  may  either  fix  the  number  of  particles  and  have  the  embedded  JavaScript  routine  calculate  the 
w/c,  or  fix  the  w/c  with  the  number  of  particles  being  calculated.  In  addition,  they  may  change 
any  of  the  specific  entries  for  the  number  of  particles  in  the  last  column  of  the  table.  For  example, 
the  number  of  particles  for  all  diameters  larger  than  41  //m  could  be  manually  set  to  zero.  For 
the  automatic  calculation  of  the  number  of  particles  when  a w/c  is  specified,  it  may  not  always  be 
possible  to  match  exactly  the  desired  w/c  ratio,  as  one  too  many  large  particles  may  result  in  a 
substantial  change  in  the  w/c.  In  this  case,  the  user  may  manually  set  the  mass  fraction  of  these 
larger  diameter  particles  to  zero  and  repeat  the  automatic  calculation  by  respecifying  the  desired 
w/c  value.  As  a final  entry,  the  user  supplies  the  name  of  the  file  to  be  created  (with  an  extension 
of  .psd  being  automatically  added  by  the  VCCTL  program). 
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Inputs  for  Program  Creatpsd 

Below  this  table,  enter  the  total  number  of  particles  or  the  desired  w/c  ratio.  You  can  also 
change  the  mass  fractions  or  specifically  fix  the  number(s)  of  certain  size  particles. 


Diameter  (pm) 

Mass  Fraction 

Number  of  particles 

1 

0 . 1348 

58217 

3 

0 . 1303 

2962 

5 

0 .09456 

504 

7 

0 . 07349 

177 

9 

0 .06085 

68 

11 

0 .0507 

30 

13 

0 . 04221 

15 

15 

0 . 03775 

9 

17 

0 . 033935 

6 

19 

0 . 028 

3 

21 

0 . 02565 

2 

23 

0 . 0233 

2 

25 

0 . 02115 

1 

27 

0.019 

1 

29 

0.019 

1 

31 

0 . 0228 

1 

35 

0 . 0347 

1 

41 

0 . 02658 

0 

47 

0 . 0454 

0 

61 

0 . 0358 

0 

73 

0 . 040025 

0 

Total  number  of  particles:  w/c  ratio  is: 

62000  .49234 

Name  to  be  given  to  this  particle  datafile:  (with  no  file  extension—  .psd  will  be  automatically 

testpart 

added) 

User  comment: 


Submit  form  to  store  created  PSD 


Return  to  the  main  menu 

Figure  5:  Form  to  specify  and  store  a PSD. 
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2.4  Menu  Selection  3)  Generate  an  initial  3-D  microstructure 

Menu  Selection  3 is  used  to  create  an  initial  3-D  microstructure  consisting  of  cement  (gypsum, 
and  filler)  particles  in  water.  For  this  form,  the  user  will  find  several  entries  that  are  highlighted 
(in  blue).  Clicking  on  any  of  these  will  direct  the  user  to  supplementary  help  information  screens. 
Often,  the  articles  referenced  in  these  help  screens  are  directly  accessible  online  as  part  of  the 
Electronic  Monograph:  Modelling  and  Measuring  the  Structure  and  Properties  of  Cement-Based 
Materials  (available  at  http://ciks.cbt.nist.gov/monograph)  [8]. 

The  form  for  Menu  Selection  3 is  shown  in  Figures  6 through  8.  The  user  must  provide  a 
negative  integer  random  number  seed  (in  the  range  [-1,-32767]),  the  thickness  of  the  single  aggregate 
particle  (slab)  placed  in  the  middle  of  the  3-D  microstructure  (if  any),  and  to  what  degree  (if  any) 
the  particles  are  to  be  dispersed  or  flocculated.  In  addition,  they  must  supply  information  on  the 
calcium  sulfate  (gypsum)  content  and  form  (dihydrate,  hemihydrate,  and  anhydrite).  Here,  the 
user  specifies  first  the  total  volume  fraction  of  randomly  distributed  calcium  sulfate  and  then  the 
fractions  of  this  total  that  are  to  be  assigned  as  hemihydrate  and  anhydrite,  respectively  (with  the 
remainder  being  distributed  as  the  dihydrate  form  of  calcium  sulfate).  The  specific  PSD  to  be  used 
in  creating  the  3-D  microstructure  must  next  be  specified.  For  the  PSD,  the  user  has  the  choice  of 
directly  specifying  the  number  and  phase  of  each  size  particle  via  the  table  in  Figure  7,  or  reading 
in  this  information  from  a PSD  datafile  created  using  Menu  Selection  2,  as  discussed  above.  Finally, 
the  user  needs  to  specify  a filename  in  which  to  store  the  created  microstructure  and  their  e-mail 
address.  Once  the  execution  of  the  genpartnew  program  is  complete,  an  e-mail  notification  will  be 
sent  to  the  user  at  the  specified  e-mail  address.  If  the  genpartnew  program  is  unable  to  place  the 
requested  numbers  of  particles,  the  user  will  be  notified  of  this  by  e-mail  as  well,  and  instructed  to 
modify  their  inputs  and  resubmit  the  form. 

Calcium  sulfate  (gypsum)  additions  can  be  handled  in  one  of  two  ways  during  the  microstructure 
creation  step.  If  only  the  composite  PSD  of  the  cement  and  gypsum  is  known,  a specific  fraction 
of  the  cement  particles  can  be  randomly  assigned  to  be  gypsum.  Conversely,  if  both  the  PSDs  of 
the  cement  and  gypsum  are  available  separately,  the  PSD  table  can  be  used  to  provide  the  specific 
numbers  of  each  size  particle  which  are  to  be  cement  and  gypsum,  respectively.  In  this  case,  the 
“Total  randomly  distributed  calcium  sulfate  fraction”  in  Figure  6 would  be  specified  to  be  0.0. 
The  same  specifications  can  be  done  with  cement  and  fly  ash  particles,  or  cement  and  inert  filler 
particles.  While  the  user  can  change  the  individual  particle  diameter  entries  in  the  table  shown 
in  Figure  7,  they  should  be  sure  that  the  entries  are  always  in  descending  order,  starting  with  the 
largest  diameter  particles.  This  is  to  ensure  that  during  random  particle  placement,  there  is  the 
best  chance  that  there  is  some  open  space  available  for  the  last  (smallest)  particles  to  be  placed. 
If  the  smallest  particles  are  randomly  located  first,  there  may  not  remain  any  empty  spaces  of 
adequate  size  to  fit  the  (last)  larger  ones. 

In  addition,  it  should  be  noted  that  only  particles  3 pixels  (pm)  in  diameter  and  larger  are 
normally  added  during  the  execution  of  Menu  Selection  3.  The  smaller  1-pixel  diameter  particles  of 
the  various  cement  clinker  phases,  gypsum,  silica  fume,  etc.  should  instead  be  specified  just  prior 
to  the  hydration  using  a fill-in  table  with  Menu  Selection  13  (see  Figure  15).  This  is  done  to  allow 
the  optional  flocculation  of  the  3 pixel  diameter  and  larger  particles  during  the  creation  of  the  3-D 
particle  microstructure  image.  Because  the  1-pixel  particles  are  generally  quite  large  in  number, 
their  inclusion  in  the  flocculation  algorithm  would  result  in  dramatic  increases  in  the  memory  and 
time  requirements  of  the  3-D  microstructure  creation  program. 
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Inputs  for  Program  Genpartnew 


If  an  entry  is  highlighted  in  blue,  you  may  click  on  it  to  receive  help  and  reference 

information! 


Random  number  seed  (negative  integer): 

-3034 


Aggregate  present  (default  is  off)  Flocculation  (default  is  none) 


Aggregate  thickness  (pixels):  Number  of  floes: 

2 1 


Dispersion  distance  between  particles  during  placement: 

0 


Sulfate  additions: 


Total  randomly  distributed  calcium  sulfate  fraction  is  with: 

0 . 05 


Hemihydrate  and  Anhydrite  fractions (on  a total  sulfate  basis; 

0.00  0.00 


Particle  Size  Distribution: 


Figure  6:  Form  to  create  a starting  3-D  particle  microstructure. 
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Input  PSD  in  table  provided  below 


Read  PSD  from  datafde  (include  extension  in  name): 

par tOO 1 . psd 


Click  on  boxes  under  Phase  to  select 


(defaults  for  cement  will  give  a w/c  of  about  0.4) 


Diameter  (pm) 

Number  of  Particles 

Phase 

41 

0 

cement(C3S) 

41 

0 

fly  ash 

37 

0 

cement(C3S) 

37 

0 

fly  ash 

35 

1 

cement(C3S) 

35 

0 

fly  ash 

33 

0 

cement(C3S) 

33  | 

0 

fly  ash 

31 

1 

cement(C3S) 

3 1 

0 

fly  ash 

29 

1 

cement(C3S) 

29 

0 

fly  ash 

27 

1 

cement(C3S) 

27 

0 

fly  ash 

25  | 

1 

cement(C3S) 

25 

0 

fly  ash 

23  | 

2 

cement(C3S) 

23  l 

0 

fly  ash 

21  l 

3 

cement(C3S) 

21 

0 

fly  ash 

19 

4 

cement(C3S) 

19  1 

0 

fly  ash 

17 

7 

cement(C3S) 

17  1 

0 

fly  ash 

15 

10 

cement(C3S) 

15  I 0 

fly  ash 

13 

18 

cement(C3S) 

13 

0 

fly  ash 

Figure  7:  Form  to  create  a starting  3-D  particle  microstructure  (continued). 
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11 

34 

cement(C3S) 

11 

0 

fly  ash 

9 

77 

cement(C3S) 

9 

0 

fly  ash 

7 

203 

cement(C3S) 

7 

0 

fly  ash 

5 

578 

cement(C3S) 

5 

0 

fly  ash 

3 

6903 

cement(C3S) 

3 

0 

fly  ash 

Name  to  be  given  to  this  microstructure:  (do  not  add  a fde  extension  as  .img  will  be  automatically 

newcem 

added) 

E-mail  address: 


Please  be  patient,  may  take  3-5  minutes  to  execute,  once  you  press  submit. 
You  will  be  e-mailed  at  the  above  address  when  execution  is  complete. 

Submit  form  to  execute  genpartnew 


Return  to  the  main  menu 


Figure  8:  Form  to  create  a starting  3-D  particle  microstructure  (continued). 
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2.5  Menu  Selection  4)  View  a 2-D  slice  from  any  3-D  microstructure 


Menu  Selection  4 is  used  to  create  and  view  a 2-D  slice  from  any  existing  3-D  microstructure 
(starting,  hydrated,  leached,  etc.).  The  user  simply  supplies  the  name  of  the  image  to  be  displayed, 
the  z-plane  (slice)  to  be  used  in  generating  the  2-D  x-y  slice  from  the  3-D  microstructure,  and  a 
magnification  factor  (IX  to  10X)  for  viewing  the  2-D  slice.  The  VCCTL  will  locate  the  proper 
slice  of  the  selected  image  and  display  a 2-D  color  image  of  the  microstructure,  along  with  a key 
identifying  the  phases  corresponding  to  each  color.  Example  images  are  provided  later  in  this 
manual,  in  Figures  23  to  26.  If  this  selection  is  used  repeatedly  with  the  same  3-D  image  filename, 
depending  on  the  user's  Internet  browser  settings  and  preferences,  they  may  need  to  use  the  Reload 
or  Refresh  button  on  their  browser  to  refresh  the  displayed  image  to  be  the  currently  selected  one. 
In  the  extreme  cases  where  this  button  doesn’t  load  the  newly  generated  image,  the  user  will  need 
to  restart  their  Internet  browser  or  use  some  other  measure  to  force  it  to  empty  the  previous  image 
from  its  cache. 

2.6  Menu  Selection  5)  Generate  an  animation  from  any  3-D  microstructure 

An  alternate  method  of  viewing  a microstructure  is  provided  by  Menu  Selection  5.  This  selection 
can  be  used  to  create  an  animated  gif  (“fly  through")  file  for  a user-specified  3-D  microstructure. 
For  this  selection,  the  user  supplies  the  following  information  before  submitting  the  displayed  form: 

1)  the  name  of  the  microstructure  to  be  animated, 

2)  a new  filename  in  which  to  store  the  animated  gif  file, 

3)  the  number  of  2-D  slices  to  include  in  the  animation  (0,100),  and 

4)  the  magnification  (1X-10X)  to  use  for  the  animation  images. 

The  first  slice  of  the  microstructure  is  always  used  as  the  starting  slice  for  the  3-D  animation 
generation.  Typically,  to  generate  a complete  (100  slice)  animation  of  a 3-D  microstructure  will 
require  15  min  to  20  min  of  CPU  time.  Thus,  the  user  should  wait  for  this  period  of  time  before 
trying  to  view  the  animation  using  Menu  Selection  C.  The  VCCTL  is  set  up  to  create  only 
one  animation  at  any  specific  time.  Therefore,  if  the  system  is  busy  creating  one 
animation  when  another  request  is  submitted,  the  user  will  be  notified  to  resubmit 
their  job  at  a later  time. 

2.7  Menu  Selection  6)  View  an  animation  from  any  3-D  microstructure 

Menu  Selection  6 is  simply  used  to  view  the  animated  gif  files  created  using  Menu  Selection  5.  The 
user  need  only  supply  the  name  of  the  animated  gif  file  to  be  displayed.  The  animated  gif  file  will 
be  displayed  five  times  from  start  to  finish,  showing  each  slice  processed  in  creating  the  animation. 

2.8  Menu  Selection  7)  Distribute  all  phases  amongst  cement  particles 

Once  an  initial  arrangement  of  cement  (and  gypsum,  fly  ash,  etc.)  particles  in  water  has  been 
created,  Menu  Selection  7 may  be  used  to  distribute  the  four  major  cement  clinker  phases  amongst 
the  cement  particles.  To  do  this,  the  3-D  microstructure  image  is  filtered  using  the  correlation  files 
measured  on  the  actual  2-D  SEM  images  of  the  cement  of  interest  [1,  2,  6].  The  filtering  process 
and  subsequent  thresholding  and  “sintering”  are  conducted  in  such  a way  as  to  match  the  volume 
and  surface  area  fractions  previously  determined  for  each  of  the  four  major  cement  clinker  phases, 
and  provided  as  input  via  the  table  shown  in  the  form  in  Figure  9.  For  example,  these  3-D  phase 
volume  and  surface  area  fractions  might  be  obtained  from  the  2-D  area  and  perimeter  fractions 
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given  in  the  corresponding  About  this  cement  page  in  the  cement  images  database  (Menu  Selection 
1).  The  user  could  simply  print  this  page  and  type  the  appropriate  values  into  the  table  for  volume 
and  surface  area  fractions. 

To  use  Menu  Selection  7,  the  user  must  specify  a negative  integer  to  be  used  as  a random 
number  seed  and  the  names  of  the  original  and  final  (processed)  microstructures,  and  select  the 
correlation  files  (from  a pulldown  menu)  to  be  used  in  filtering  the  original  microstructure.  In 
addition,  they  must  fill  in  the  table  providing  the  volume  and  surface  area  fractions  of  each  of  the 
four  major  cement  clinker  phases  and  provide  their  e-mail  address  in  order  to  receive  a notification 
when  the  processing  is  complete.  Ordinarily,  on  the  order  of  20  min  of  CPU  time  is  required  to 
distribute  the  various  phases  amongst  the  original  cement  particles.  The  VCCTL  is  set  up  to 
perform  only  one  phase  distribution  at  any  specific  time.  Therefore,  if  the  system  is 
busy  distributing  phases  for  one  microstructure  when  another  request  is  submitted, 
the  user  will  be  notified  to  resubmit  their  job  at  a later  time. 
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Inputs  for  Program  Distrib3d 


Access  the  latest  version  of  the  hydration  manual  to  read  about  phase  distribution 


Random  number  seed  (negative  integer): 

-3034 


Filename  of  original  microstructure: 

cem!33wc030nl . img 


Filename  of  the  final  microstructure: 

ceml33wc030nlf . img 

Processed  cement  to  use  for  the  correlation  files: 

ceml33r 

Phase  volume  (area)  and  surface  fractions  for  the  cement  being  generated: 


Phase 

Volume  fraction 

Surface  area 
fraction 

c3s 

0.7018 

0 . 6491 

c2s 

0 . 1315 

0 . 1764 

c3a 

0 . 0827 

0 . 1138 

c4af 

0 . 0840 

0 .0607 

E-mail  address: 

You  will  be  e-mailed  at  the  above  address  when  execution  is  complete. 

Submit  form  to  execute  distrib3d 

Return  to  the  main  menu 

Figure  9:  Form  to  specify  inputs  to  be  used  in  distributing  phases  amongst  cement  particles. 
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2.9  Menu  Selection  8)  Filter  a 3-D  microstructure  to  distribute  only  two  phases 
amongst  cement  particles 

Normally,  the  user  will  apply  Menu  Selection  7 to  automatically  distribute  all  four  phases  amongst 
the  initial  cement  particles.  However,  in  those  cases  where  special  control  over  the  distribution 
process  is  needed,  the  user  may  control  each  step  in  the  distribution  process  using  Menu  Selections 
8,  9,  and  10.  Specifically,  Menu  Selection  8 is  used  to  filter  a 3-D  image  using  a specific  correlation 
file  for  a specific  cement.  Basically,  an  image  of  random  noise  is  superposed  on  the  random  particle 
microstructure  image  and  filtered  by  (convolved  with)  the  correlation  function  for  one  or  more 
phases  determined  from  the  2-D  SEM  images.  Then,  the  correlated  noise  image  is  thresholded 
to  separate  one  phase  of  the  microstructure  into  two  phases,  with  the  user  specifying  the  desired 
volume  fraction  of  each  phase.  For  example,  initially,  the  cement  particles  are  separated  into  silicates 
and  aluminates  using  a cemXXX.sil  correlation  file  [1].  In  subsequent  passes  of  this  algorithm,  the 
silicates  are  separated  into  C3S  and  C2S , while  the  aluminates  are  separated  into  C3A  and  C4AF. 
When  the  user  executes  Menu  Selection  7,  all  three  passes  are  conducted  in  one 
execution  of  the  program. 

To  execute  the  filtering  ( rand3d ) program,  the  user  must  supply  a random  number  seed,  the 
names  of  the  input  and  output  microstructures,  the  correlation  file  to  be  used,  the  volumetric  phase 
fraction  of  the  first  phase  to  be  retained  after  the  filtering,  and  their  e-mail  address.  As  with  many 
of  the  other  menu  selections,  the  user  will  be  notified  by  an  e-mail  message  once  execution  of  the 
filtering  is  complete.  The  VCCTL  is  set  up  to  perform  only  one  microstructure  filtering 
at  any  specific  time.  Therefore,  if  the  system  is  busy  filtering  one  microstructure  when 
another  request  is  submitted,  the  user  will  be  notified  to  resubmit  their  job  at  a later 
time. 


16 


Inputs  for  Program  Rand3d 


Access  the  latest  version  of  the  hydration  manual  to  read  about  image  filtering 


Random  number  seed  (negative  integer): 

-3034 


Filename  of  original  microstructure: 

newcem. img 


Filename  of  the  new  microstructure: 

newcemla . img 

Filename  of  the  correlation  file: 

ceml33r.si! 

Fraction  of  first  phase  to  be  maintained  as  is: 

0.50 


E-mail  address: 


You  will  be  e-mailed  at  the  above  address  when  execution  is  complete. 

Submit  form  to  execute  rand3d 

Return  to  the  main  menu 

Figure  10:  Form  to  specify  inputs  to  be  used  in  filtering  a 3-D  microstructure. 
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2.10  Menu  Selection  9)  Compute  spatial  statistics  for  a 3-D  microstructure 

Menu  Selection  9 can  be  used  to  measure  the  volume  and  surface  pixel  counts  for  each  phase 
present  in  any  3-D  microstructure.  This  menu  selection  can  be  employed  during  the  creation  of  an 
initial  microstructure,  after  hydration  of  an  initial  microst.ructure,  or  after  leaching  of  a hydrated 
microstructure.  To  use  this  selection,  the  user  need  only  supply  the  name  of  the  3-D  microstructure 
image  file  to  be  analyzed.  Within  one  minute  to  two  minutes  from  submitting  the  form,  a table  of 
the  volume  and  surface  counts,  such  as  that  shown  in  Figure  11,  will  be  returned  to  the  user.  In  this 
case,  the  volume  and  surface  fractions  for  the  four  major  clinker  phases  in  an  analyzed  image  can 
be  compared  to  those  requested  using  Menu  Selection  7 (see  Figure  9),  with  a very  good  agreement 
being  observed. 
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Resultant  statistics  for  cenil33wc030nlfimg 


Phase 

Phase 

Volume 

Surface 

Volume 

Surface 

ID 

count 

count 

fraction 

fraction 

Porosity 

0 

569660 

0 

C3S 

1 

284464 

293128 

0 . 69948 

0 . 64804 

C2S 

2 

53645 

80013 

0 . 13191 

0 . 17689 

C3A 

3 

33982 

51716 

0.08356 

0 .11433 

C4AF 

4 

34591 

27471 

0.08506 

0 .06073 

Total 

406682 

452328 

Gypsum 

5 

23658 

23382 

Hemihyd . 

6 

0 

0 

Anhydrite 

7 

0 

0 

Pozzolan 

8 

0 

0 

Inert 

9 

0 

0 

ASG 

10 

0 

0 

CAS2 

11 

0 

0 

CH 

12 

0 

0 

C-S-H 

13 

0 

0 

C3AH6 

14 

0 

0 

Ettringite 

15 

0 

0 

Ettringi te 

16 

0 

0 

Afm 

17 

0 

0 

FH3 

18 

0 

0 

PozzC-S-H 

19 

0 

0 

CaC12 

20 

0 

0 

Freidelsal t 

21 

0 

0 

Strati ingite2 2 

0 

0 

Sec . gypsum 

23 

0 

0 

Inertagg 

24 

0 

0 

Abs . gypsum 

25 

0 

0 

Fly  ash 

26 

0 

0 

C3A ( fly  ash) 3 5 

0 

0 

Return  to  the  main  menu 


Figure  11:  Results  from  spatial  statistics  analysis  of  an  initial  3-D  microstructure. 
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2.11  Menu  Selection  10)  Sinter  a 3-D  microstructure  to  adjust  the  hydraulic 
radius  of  one  phase 


Following  the  filtering  process  (Menu  Selection  8),  the  surface  area  fractions  of  the  filtered  phases 
(from  Menu  Selection  9)  may  not  match  those  of  the  phases  in  the  starting  2-D  SEM  images.  These 
surface  area  fractions  can  be  adjusted  using  a 3-D  “sintering”  algorithm  to  swap  specific  locations 
of  the  two  phases  in  a manner  that  obtains  the  desired  surface  area  fractions.  Details  of  the  3- 
D sintering  algorithm  have  been  published  previously  [9],  with  the  specific  application  to  cement 
particles  described  in  the  CEMHYD3D  Version  2.0  user’s  manual  [1]. 

The  input  form  for  Menu  Selection  10  is  shown  in  Figure  12.  The  user  must  supply  a negative 
integer  random  number  seed,  the  names  of  the  input  and  output  images,  the  desired  (target) 
hydraulic  radius  of  the  phase  being  sintered,  and  the  phase  IDs  of  the  two  phases  between  which 
to  execute  the  sintering  algorithm.  The  target  hydraulic  radius  of  a phase  can  be  calculated  (based 
on  the  volume  count  provided  with  Menu  Selection  9)  as  [1]: 


volume  of  phase  in  pixels 
desired  surface  area  count  in  pixels 


(1) 


The  factor  of  6 over  4 is  included  to  correct  for  the  difference  between  the  approximate  surface  area 
of  a digitized  (pixel-based)  sphere  and  a continuum  one  of  the  same  diameter.  It  should  be  noted 
that  the  hydraulic  radius  of  a phase  can  only  be  increased  when  applying  the  sintering  algorithm  to 
it.  If  it  is  necessary  to  decrease  the  hydraulic  radius  of  a phase,  the  algorithm  should  be  executed 
so  as  to  increase  the  hydraulic  radius  of  the  secondary  phase,  thus  decreasing  that  of  the  primary 
one  [1]. 
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Sinter  a 3-D  microstructure  to  adjust  hydraulic 

radius 


Please  refer  to  the  cement  hydration  manual  for  more  details. 


Random  number  seed: 

-24 


Image  file  name  to  be  sintered(complete  with  extension): 

input . img 


Desired  hydraulic  radius  for  first  phase: 

1.50 


Phases  to  execute  sintering  between:  and 

14 


New  image  file  name  to  create(complete  with  extension): 

output . img 


Please  be  patient,  may  take  3-5  minutes  to  sinter  the  image  after  you  submit. 

Submit  form  to  sinter  microstructure 

Return  to  the  main  menu 

Figure  12:  Form  to  specify  inputs  to  be  used  in  “sintering”  a 3-D  microstructure. 
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2.12  Menu  Selection  11)  Distribute  fly  ash  phases  amongst  fly  ash  particles 

While  Menu  Selections  8 to  10  could  possibly  be  used  to  distribute  phases  amongst  fly  ash  particles 
as  well  as  cement  particles,  Menu  Selection  1 1 has  been  developed  specifically  for  the  distribution 
of  phases  amongst,  fly  ash  particles,  under  two  simplifying  cases.  In  the  first  case,  the  particles  are 
assumed  to  be  monophase  and  the  particles  are  randomly  assigned  phases  to  achieve  the  desired 
phase  volume  fractions.  In  the  second,  the  phases  are  simply  randomly  distributed  amongst  all  of 
the  pixels  identified  as  fly  ash  in  the  3-D  particle  microstructure.  This  would  imply  a very  fine  and 
relatively  uniform  distribution  of  the  fly  ash  phases  amongst  the  individual  particles.  Currently, 
the  following  phases  are  considered  for  distribution  amongst  the  fly  ash  particles:  aluminosilicate, 
calcium  aluminodisilicate,  silica,  anhydrite,  calcium  chloride,  tricalcium  aluminate,  and  an  inert 
phase.  Reactions  between  these  phases  and  the  main  components  of  cement  have  been  incorporated 
into  the  current  version  of  the  cement  hydration  model  [1,  10]. 

To  execute  Menu  Selection  11,  the  user  provides  a negative  integer  random  number  seed,  selects 
whether  distribution  is  to  be  done  on  a random  particle  or  a random  pixel  basis,  and  fills  in  the 
remainder  of  the  form  shown  in  Figure  13.  The  user  must  provide  the  original  microstructure  and 
particle  microstructure  filenames,  the  name  of  an  output  file  to  hold  the  processed  microstructure, 
the  number  of  fly  ash  pixels  in  the  3-D  microstructure  (can  be  determined  using  Menu  Selection  9 
above),  the  desired  volume  fractions  for  the  various  fly  ash  phases,  and  their  e-mail  address.  The 
particle  microstructure  image  is  needed  for  the  case  where  the  phases  are  distributed  randomly 
amongst  the  particles,  one  phase  per  particle.  When  fly  ash  particles  are  placed  during  the  creation 
of  a microstructure  (Menu  Selection  3),  they  are  assigned  a phase  ID  of  26.  Thus,  when  using 
Menu  Selection  9,  the  user  must  be  sure  to  note  the  volume  count  for  phase  ID  26  and  use  it  as 
the  correct  value  for  the  fly  ash  count  in  the  form  shown  in  Figure  13. 
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Inputs  for  Distributing  Phases  within  Fly  Ash 


Random  number  seed  (negative  integer):- 32 17 
Distribute  fly  phases  randomly  on  a particle  basis 
Distribute  fly  phases  randomly  on  a pixel  basis 
Name  of  input  phase  microstructure  file  (with  extension) : 

newcem. img 

Name  of  input  particle  microstructure  file  (with  extension) : 

pnewcem. img 

(typically  this  is  the  phase  microstructure  file  with  a p in  front  of  it) 


Name  of  output  phase  microstructure  file  to  create  (with  extension) : 

newceml . img 

Number  of  fly  ash  pixels  (from  executing  stat3d-  phase  26  is  fly  ash) 

10000 

Phase  fractions  (remainder  will  be  considered  as  inert): 

Phase  Volume  (number)  fraction 

Aluminosilicate  glass: 


Calcium  alumina  disilicate 

Tricalcium  aluminate: 

Calcium  chloride: 

Silica : 

Anhydrite  (CaS04): 


0 . 10 

0 .25 

0 . 00 

0 . 00 
0 . 40 

0 . 05 


E-mail  address: 

dale . bent z@nist . gov 


Submit  form  to  execute  distribution  of  fly  ash  phases 

Return  to  the  main  menu 

Figure  13:  Form  to  specify  inputs  to  be  used  in  distributing  the  phases  amongst  the  fly  ash  particles. 
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2.13  Menu  Selection  12)  Input  a temperature  history  (curing)  file 

The  CEMHYD3D  cement  hydration  model  can  be  executed  under  a variety  of  curing  conditions 
(isothermal,  adiabatic,  or  user- programmed  temperature  history).  Menu  Selection  12  is  used  to 
specify  a controlled  temperature  history  to  be  used  in  the  hydration  process  (Menu  Selection  13). 
The  user  simply  fills  in  a table  of  starting  and  ending  times  with  the  corresponding  temperatures. 
A linear  change  in  temperature  with  time  is  assumed  between  the  starting  and  ending  points. 
Currently,  the  table  is  limited  to  six  separate  entries  which  should  be  adequate  to  specify  a steam 
curing  profile,  for  example.  In  addition  to  the  values  for  the  table,  the  user  also  specifies  a filename 
in  which  to  store  the  temperature  history  so  that  it  may  be  accessed  later  during  the  execution  of 
the  hydration  module. 

Inputs  for  Program  Tempstore 

Filename  to  be  created: 

tempnew . dat 

Temperature  History  for  Curing: 

Note  that  the  starting  and  ending  times  should  define  a piecewise  continuous  curve  (e.g.,  starting  time 

(1+1)  should  equal  ending  time  (I)). 


Starting  time 

(h) 

Ending  time  (h) 

Starting  temperature  (C) 

Ending  temperature 

(C) 

0 . 0 

4 . 0 

25.0 

25.0 

4.0 

5.0 

25.0 

60.0 

5 . 0 

29 . 0 

60.0 

60.0 

29.0 

40 . 0 

60.0 

25.0 

40 . 0 

1000 . 0 

25.0 

25.0 

1000.0 

40000 . 0 

25.0 

25.0 

Submit  form  to  store  temperature  history 


Return  to  the  main  menu 

Figure  14:  Form  to  specify  the  temperature  history  to  be  used  in  curing  a microstructure  during 
hydration. 
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2.14  Menu  Selection  13)  Hydrate  a 3-D  microstructure 

The  core  of  the  VCCTL  is  Menu  Selection  13  as  it  permits  one  to  read  in  a starting  microstructure, 
add  some  one-pixel  particles  (if  needed),  and  perform  hydration  under  user-specified  curing  condi- 
tions. The  fill-in  forms  for  doing  this  are  shown  in  Figures  15  through  17.  As  with  several  other 
selections  from  the  main  menu,  the  VCCTL  is  currently  set  up  to  perform  only  one 
hydration  at  any  specific  time.  Therefore,  if  the  system  is  busy  executing  a hydration 
when  another  request  is  submitted,  the  user  will  be  notified  to  resubmit  their  job  at 
a later  time. 

To  begin  with,  the  user  supplies  a negative  integer  random  number  seed,  and  the  names  of  the 
starting  cement  microstructure  and  particle  image  (typically  the  original  image  filename  with  the 
letter  p in  front  of  it)  files  to  be  used  in  the  hydration.  Then,  they  have  the  option  of  adding  some 
one-pixel  particles  of  various  phases  prior  to  the  execution  of  the  hydration.  This  table  should 
be  used  to  add  the  one-pixel  (1  pm  diameter)  particles  necessary  to  complete  the  cement  PSD 
[1].  The  exact  number  of  one-pixel  particles  for  each  of  the  four  cement  phases  and  gypsum  (or 
hemihydrate  or  anhydrite)  can  be  determined  from  the  volume  counts  (Menu  Selection  9)  for  the 
final  created  3-D  microstructure  along  with  the  PSD  established  via  Menu  Selection  2,  for  example. 
Menu  selection  2 will  provide  the  total  number  of  one-pixel  particles  to  add,  while  Menu  Selection 
9 (along  with  the  known  phase  compositon  of  the  cement)  will  allow  calculation  of  how  this  total 
number  of  one-pixel  particles  should  be  distributed  to  create  a 3-D  system  that  matches  the  desired 
overall  phase  composition  of  the  cement  [1].  A simple  JavaScript  routine  to  aid  in  this  calculation 
is  available  as  a help  screen  by  clicking  on  the  highlighted  Number  of  one  pixel  particles  to  add:. 
Additionally,  one-pixel  particles  of  silica  fume,  an  inert  filler,  or  one  phase  of  the  user’s  choice  can 
be  added  via  this  table. 

The  user  then  specifies  the  parameters  that  will  control  the  hydration  including: 

1)  the  number  of  cycles  of  hydration  to  be  executed, 

2)  the  thermal  curing  conditions  (isothermal,  adiabatic,  or  temperature-defined), 

3)  the  initial  saturation  curing  conditions  (saturated  or  sealed), 

4)  whether  the  conversion  of  primary  C—S—H  to  pozzolanic  C—S—H  is  prohibited  or  allowed 
(this  selection  is  only  relevant  when  pozzolanic  materials  such  as  silica  fume  or  fly  ash  are 
present  in  the  starting  microstructure), 

5)  the  frequencies  in  cycles  to  evaluate  the  water-filled  porosity  percolation,  the  percolation 
of  total  solids,  and  the  individual  particle  hydration  rates  (these  can  be  “disabled”  by  setting 
the  frequencies  to  values  larger  than  the  requested  number  of  cycles  of  hydration), 

6)  the  initial  curing  temperature, 

7)  an  estimate  of  the  activation  energy  for  cement  hydration, 

8)  an  estimate  of  the  activation  energy  for  the  pozzolanic  reactions, 

9)  the  conversion  factor  to  go  from  cycles  of  hydration  to  h of  real  time,  and 

10)  the  aggregate  volume  fraction  in  a concrete  mixture  (used  for  prediction  of  the  adiabatic 
heat  signature  of  a concrete  — can  be  set  to  0.0  or  left  as  is  for  cement  paste  systems). 

Finally,  the  user  once  again  supplies  their  e-mail  address  so  that  they  may  be  notified  when  the 
hydration  is  completed.  All  of  the  parameters  listed  above  are  described  in  detail  in  Version  2.0 
of  the  CEMHYD3D  user’s  manual  [1],  and  some  of  the  entries  in  the  form  are  highlighted  in 
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blue,  indicating  that  supplemental  help  screens  are  directly  available  for  viewing  by  clicking  on  the 
highlighted  blue  entry. 

During  the  hydration,  a number  of  output  fdes  are  created  by  the  CEMHYD3D  program.  These 
all  utilize  the  starting  microstructure  filename  as  their  root  and  add  on  appropriate  extensions  to 
indicate  the  results  that  each  one  contains.  For  example,  if  a file  cemnew.img  is  hydrated  for 
4000  cycles  of  hydration  at  a starting  temperature  of  25  °C  and  under  isothermal  and  saturated 
conditions,  the  following  output  files  will  normally  be  created: 

1)  cemnew.img. 4000. 25. 000  - microstructure  image  file  after  the  hydration;  the  4000  indicates 
the  number  of  cycles  of  hydration,  the  25  indicates  the  starting  temperature,  and  the  000 
indicates  the  selections  for  conversion  of  primary  C—S—H  to  pozzolanic  C—S—H  (prohibited^, 
allowed=l),  thermal  curing  (isothermal^O,  adiabatic=l,  temperature-defined=2),  and  initial 
saturation  (saturated=0,  sealed=l), 

2)  cemnew. heat. 4000. 25. 000  - output  file  containing  a listing  of  the  degree  of  hydration  (both 
volume  and  mass  bases),  computed  heat  released,  and  computed  gel-space  ratio  [11]  vs.  num- 
ber of  hydration  cycles  and  computed  real  hydration  time, 

3)  cemnew. adi. 4000. 25. 000  - output  file  containing  a listing  of  the  system  temperature  vs. 
computed  real  hydration  time, 

4)  cemnew. chs. 4000. 25. 000  - output  file  containing  a listing  of  the  chemical  shrinkage  vs. 
number  of  hydration  cycles  and  computed  real  hydration  time, 

5)  cemnew. pha. 4000. 25. 000  - output  file  containing  a listing  of  the  phase  volume  counts  vs. 
number  of  hydration  cycles, 

6)  cemnew. pps. 4000. 25. 000  - output  file  containing  a listing  of  the  results  for  the  examination 
of  the  percolation  of  the  water-filled  capillary  porosity  vs.  cycles,  computed  hydration  time, 
and  degree  of  hydration, 

7)  cemnew. pts. 4000. 25. 000  - output  file  containing  a listing  of  the  results  for  the  examination 
of  the  percolation  of  the  total  solids  vs.  cycles,  computed  hydration  time,  and  degree  of 
hydration,  and 

8)  cemnew. phr. 4000. 25. 000  - output  file  containing  a listing  of  the  results  for  the  hydration 
degree  of  individual  particles  after  specific  numbers  of  hydration  cycles. 

When  Menu  Selection  13  is  executed,  all  of  the  user  input  parameters  are  echoed  back  in  a form, 
along  with  the  names  of  all  of  the  output  files  that  may  be  created  by  the  CEMHYD3D  program. 
However,  normally,  users  don’t  need  to  concern  themselves  with  the  contents  of  these 
files,  as  they  can  instead  use  Menu  Selection  17  to  plot  any  specific  property  of  interest 
for  the  systems  that  have  been  hydrated  using  Menu  Selection  13. 
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Inputs  for  Program  CEMHYD3D 


Access  the  latest  version  of  the  hydration  manual 

Random  number  seed  (negative  integer): 

-3999 

Filename  of  starting  cement  microstructure: 

ceml3  3wc03  Onl f . img 

Filename  of  original  particle  microstructure: 

pceml3 3wc03 Onl . img 

Number  of  one  pixel  particles  to  add: 

Phase  Number  to  add 

C7S : 

45921 

C9S : 

. 17065 

Figure  15:  Form  to  hydrate  a starting  3-D  microstructure. 
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5480 


C3A: 

C4AF: 

Gypsum: 
Hemihydrate : 
Anhydrite : 
Silica  fume: 

Inert  filler: 

Phase  : 


4706 

5201 

0 

0 

0 

0 


0 0 

Number  of  cycles  of  hydration  to  execute: 


100 


Curing  thermal  conditions:  Temperature  history  filename: 

isothermal  tempone.dat 

Curing  saturation  conditions: 

saturated 

Conversion  of  primary  C-S-H  to  pozzolanic  C-S-H  is: 

prohibited 

Frequency  (in  cycles)  to  evaluate  percolation  of  porosity: 


Frequency  (in  cycles)  to  evaluate  percolation  of  total  solids  (set): 


Frequency  (in  cycles)  to  evaluate  individual  particle  hydration: 

5000 


Initial  temperature:  °C 

25 . 0 


Figure  16:  Form  to  hydrate  a starting  3-D  microstructure  (continued). 
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Activation  energy  for  cement  hydration:  kJ/mole 

40.0 


Activation  energy  for  pozzolanic  reactions:  kJ/mole 

83.14 


Conversion  factor  to  go  from  cycles  to  time: 

0.0003 


Aggregate  volume  fraction  in  a concrete  mixture: 

0 . 72 


E-mail  address: 

dale . bent z@nis t . gov 


Submit  form  to  execute  cemhyd3d 


Return  to  the  main  menu 


Figure  17:  Form  to  hydrate  a starting  3-D  microstructure  (continued). 


29 


2.15  Menu  Selection  14)  Measure  phase  fractions  vs.  distance  from  aggregate 

Menu  Selection  14  can  be  used  to  analyze  the  phase  fractions  present  as  a function  of  distance 
from  an  aggregate  surface,  for  microstructures  that  were  originally  created  containing  an  aggregate 
particle.  The  user  supplies  the  name  of  the  microstructure  to  be  analyzed,  and  a file  with  an 
.agg  extension  replacing  the  .img  portion  of  the  submitted  filename  is  created  containing  the 
quantitative  phase  count  results.  In  addition,  if  the  user  provides  their  e-mail  address,  this  same 
file  will  be  sent  to  them  for  further  analysis  and  plotting,  in  their  own  spreadsheet  package  for 
example.  For  each  incremental  pixel  distance  away  from  the  aggregate,  the  file  lists  the  distance  in 
pixels  and  the  phase  counts  for  phase  IDs  0 (water-filled  porosity)  to  26  (fly  ash)  and  phase  IDs  35 
(C3/I  in  fly  ash  before  hydration  is  executed)  and  38  (empty  porosity).  All  of  the  phases  and  their 
corresponding  IDs  are  provided  in  Figure  11. 

2.16  Menu  Selection  15)  Leach  CH/cement  from  a hydrated  microstructure 

Menu  Selection  15  can  be  used  to  selectively  leach  (remove)  a portion  or  the  complete  volume 
of  a phase(s)  from  a hydrated  microstructure.  The  calcium  hydroxide  (CH)  and  the  initial  four 
cement  clinker  phases  are  all  leachable,  with  the  user  selecting  which  phases  to  include  in  a specific 
execution  of  the  leaching  algorithm.  As  shown  in  the  input  form  in  Figure  18,  the  user  must  supply 
the  filename  of  the  microstructure  to  be  leached,  a new  filename  in  which  to  store  the  microstructure 
after  leaching,  the  phases  to  be  leached,  the  number  of  leaching  cycles  to  execute  or  to  entirely 
leach  the  selected  phases  from  the  microstructure,  a leaching  probability  (pi  in  the  range  0.0- 1.0), 
a negative  integer  random  number  seed  for  the  leaching  process,  and  their  e-mail  address.  The 
leaching  probability  can  be  used  to  control  the  amount  of  a phase(s)  leached  in  a given  cycle.  If  a 
given  leachable  phase  pixel  is  in  contact  with  pore  space,  it  will  be  leached  with  probability  pi  (and 
remain  as  is  with  probability  1-p/).  In  this  way,  microstructures  with  “leached”  levels  less  than  that 
removed  during  a single  cycle  of  full  (p/  = 1.0)  leaching  can  be  created  for  subsequent  evaluation. 
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Leaching  of  phases  from  a 3-D  microstructure 


Access  a paper  about  the  effects  of  leaching  on  pore  space  percolation  and  diffusi vity. 

Image  file  name  to  be  leached  (complete  with  extension): 

ceml33wc030nlf . img 

File  name  for  final  leached  image  (complete  with  extension): 

ceml33wc030nlf lch. img 

Phases  to  be  leached  (default  is  CH  only): 

CH  c3s  c2s  c3a  c4af 

Leaching  execution: 

Leach  for  cycles 
100 

Leach  selected  phases  completely  from  microstructure 
1 • 0 leaching  probability 
~34  random  number  seed  for  leaching 
E-mail  address: 


It  may  take  15-20  minutes  to  leach  the  image  after  you  submit. 

You  will  be  notified  by  e-mail  at  the  above  address  when  leaching  is  complete. 

Submit  form  to  execute  the  leaching 

Return  to  the  main  menu 

Figure  18:  Form  to  specify  inputs  for  leaching  CH/cement  from  a hydrated  3-D  microstructure. 
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2.17  Menu  Selection  16)  Evaluate  percolation  properties  of  a hydrated/leached 
microstructure 

Menu  Selection  16  can  be  used  to  evaluate  the  percolation  properties  of  a phase(s)  in  a hy- 
drated/leached microstructure.  To  use  this  selection,  the  user  supplies  the  filename  of  the  mi- 
crostructure to  be  evaluated  and  selects  the  phase(s)  to  be  evaluated  from  a pull-down  menu.  In 
version  1.0  of  the  VCCTL,  the  percolation  properties  of  the  following  phases  can  be  evaluated: 
water-filled  porosity,  total  porosity,  CH,  C—S—H , pozzolanic  C—S—H , total  C—S—H , ettringite,  or 
stratlingite.  An  example  of  the  results  for  a percolation  evaluation  of  the  water-filled  porosity  in 
an  initial  microstructure  (ceml33wc030nlf.img)  is  provided  in  the  Table  shown  in  Figure  19.  For 
each  of  the  three  principal  directions  (x,  y,  and  z),  the  VCCTL  returns  the  total  number  of  pixels  of 
the  phase  of  interest,  the  fraction  of  these  pixels  which  are  part  of  percolated  (through)  pathways 
through  the  microstructure,  and  the  fraction  of  the  total  pixels  which  are  accessible  from  one  face 
of  the  microstructure.  For  the  example  shown  in  Figure  19,  since  no  hydration  has  yet  occurred, 
the  initial  water-filled  porosity  between  cement  particles  is  highly  percolated,  with  nearly  all  of  the 
water-filled  porosity  pixels  being  part  of  a percolated  pathway. 

Resultant  percolation  properties  for  Water-filled  porosity  in 

ceml33wc030nlf.img 


Direction 

Total 

Percolated 

Accessible 

- 

count 

fraction 

fraction 

X 

569660 

0.9994 

0.9995 

Y 

569660 

0.9994 

0.9994 

Z 

569660 

0.9994 

0.9994 

Return  to  the  main  menu 

Figure  19:  Results  for  percolation  of  capillary  porosity  in  microstructure  ceml33wc030nlf.img. 


2.18  Menu  Selection  17)  Plot  predicted  properties  (vs.  experimental  data) 

Menu  Selection  17  is  used  both  to  plot  the  different  predicted  properties  and  to  enter  experimental 
data  for  comparison  against  the  predicted  values.  Upon  making  this  selection,  the  user  is  presented 
with  a form  with  two  pull  down  menus  to  select  either  a specific  system  for  which  to  list  the  available 
plots  or  for  which  to  enter  a set  of  experimental  data,  as  shown  in  Figure  20.  Upon  submitting  this 
first  form,  a second  form  containing  a table  of  the  various  options  available  either  for  plotting  or 
data  entry  is  returned  (see  Figure  21  for  an  example  of  the  table  returned  in  the  case  of  plotting). 

With  this  table,  in  the  case  of  plotting,  the  user  can  select  to  create  all  available  plots  for  a 
given  system,  or  simply  to  create  a single  plot  of  interest.  In  all  cases,  the  user  can  elect  to  include 
the  experimental  data  (when  available)  on  the  plot.  Typically,  the  computer  model  data  is  plotted 
using  blue  squares  connected  by  a blue  line,  while  the  experimental  data  is  shown  as  discrete  red 
diamond  data  points.  The  plot  is  typically  labelled  with  the  name  of  the  datafile  from  which  it  was 
created. 
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Virtual  Cement  and  Concrete  Testing 
Laboratory 

List  data  files  for  plotting: 

aalborgwc27f 

List  data  files  for  data  entry: 

aalborgwc27f 


Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  20899 


Figure  20:  Form  for  selecting  files  to  plot  or  for  which  to  enter  experimental  data. 


In  Version  1.0  of  the  VCCTL,  the  following  plots  are  available  (when  the  appropriate  data  has 
been  created  during  a hydration  run): 

1)  adiabatic  heat  signature  plot-  plot  of  system  temperature  vs.  computed  hydration  time 
(not  very  interesting  for  hydration  under  isothermal  conditions), 

2)  chemical  shrinkage-  plot  of  chemical  shrinkage  (ml/g  cement)  vs.  computed  hydration  time 

(h), 

3)  hydration-  plot  of  degree  of  hydration  (mass  basis)  vs.  computed  hydration  time, 

4)  heat  release-  plot  of  cumulative  heat  release  (kJ/kg  cement)  vs.  computed  hydration  time, 

5)  compressive  strength-  plot  of  estimated  compressive  strength  (MPa)  vs.  computed  hydra- 
tion time  (the  user  must  provide  a pre-factor  to  be  used  in  the  calculation  of  compressive 
strength  based  on  Power’s  gel-space  ratio  theory  [11]), 

6)  capillary  porosity  percolation-  fraction  connected  capillary  porosity  vs.  either  total  porosity 
or  computed  hydration  time,  and 

7)  total  solids  percolation-  fraction  connected  “total  solids”  vs.  either  computed  hydration 
time  or  degree  of  hydration. 

If  the  user  elects  to  enter  experimental  data  corresponding  to  a specific  simulated  system,  they 
have  the  option  of  downloading  the  data  into  the  VCCTL  from  a datafile  located  on  their  own 
computer  or  to  enter  the  data  manually.  In  the  latter  case,  the  VCCTL  will  provide  a fill-in 
table  with  two  columns  corresponding  to  the  values  to  be  plotted  on  the  x and  y axes,  respectively. 
For  example,  for  the  case  of  an  adiabatic  heat  signature,  the  user  would  have  to  enter  the  time- 
temperature  results  into  this  table  and  then  simply  submit  the  form  (using  the  “Copy  data”  button 
above  the  data  entry  table)  to  store  this  datafile  for  later  access  during  plotting.  Always,  when 
experimental  data  has  been  entered,  the  user  has  the  option  of  including  it  on  or  omitting  it  from 
the  associated  simulation  results  plot. 
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Virtual  Cement  and  Concrete  Testing 
Laboratory 

Plot  Selection  Form 


Case:  cem  1 36wc04852f 


Include  data:  ® 
strength  prefactor 

l95 iMPa 

Plot 

Type 
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Pozzolanic 

CSH 

Conversion 

Heat 

Transfer 
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mm mm  j 
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25 
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1600 
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saturated 

25 

hydration 

1600 

off 

isothermal 

saturated 

25 
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1600 
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25 

Include  data:  ® 

compressive  strength 
strength 
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1600 
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25 

Building  and  Fire  Research  Laboraton ' 
National  Institute  of  Standards  and  Technology’ 

Gaithersburg,  MD  20899 


Figure  21:  Example  of  results  available  for  plotting  via  Menu  Selection  17. 
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2.19  Menu  Selection  18)  Predict  chloride  ion  diffusivity  of  a concrete  based  on 
mixture  parameters 

The  final  menu  selection  of  the  VCCTL  is  provided  to  calculate  an  estimated  chloride  ion  diffusion 
coefficient  for  a concrete  (optionally  containing  silica  fume),  based  on  knowledge  of  its  mixture 
proportions  and  the  achieved  degree  of  hydration.  The  equations  used  in  this  prediction  have 
been  presented  previously  [12]  and  are  also  incorporated  into  a prototype  Computer  Integrated 
Knowledge  System  (CIKS)  for  predicting  the  service  life  of  steel-reinforced  concrete  exposed  to 
chloride  ions  [13],  available  separately  from  the  VCCTL.  To  use  Menu  Selection  18,  the  user  supplies 
the  following  four  inputs  for  the  concrete  of  interest:  the  w/c  ratio,  the  silica  fume  content  (mass 
fraction  of  cement  basis),  the  volume  fraction  of  aggregates,  and  the  degree  of  hydration  of  the 
cement.  Recommended  ranges  of  applicability  for  each  of  these  parameters  are  provided 
on  the  fill-in  form  and  the  user  should  extrapolate  beyond  these  ranges  only  with 
extreme  caution.  Upon  submittal,  the  estimated  value  of  the  diffusion  coefficient,  D,  in  SI  units, 
is  returned  along  with  the  90  % confidence  limits  for  D and  some  instructions  on  methodologies  to 
lower  the  D value,  if  necessary. 
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Prediction  of  the  Chloride  Ion  Diffusivity  of  a 

Concrete  Based  on  Mixture  Parameters 


Please  supply  the  following  parameters  (defaults  provided) 

w/c  ratio  (water-to-cement  ratio)  (recommended  range  of  0.3  to  0.5) 

0 .40 

Silica  fume  addition  (mass  % of  cement)  (recommended  range  of  0 % to  10  %) 

5 . 00 

Volume  fraction  of  Aggregate  % (recommended  range  of  62  % to  70  %) 

Degree  of  Hydration  of  the  Cement  (recommended  range  of  0.6  to  0.9) 

Submit  form  to  estimate  diffusivity 

Reset  all  values  to  defaults 


Papers  on  the  modelling  techniques  used  to  estimate  D are 
available  online  by  clicking  on  the  titles  in  the  following 
references : 


Bentz,  D.P.,  Influence  of  Silica  Fume  on  Diffusivity  in  Cement-Based  Materials.  II.  Multi-Scale 
Modeling  of  Concrete  Diffusivity,  Cement  and  Concrete  Research , 30  (7),  1 121-1 129,  2000. 

Bentz,  D.P.,  Detwiler,  R.J.,  Garboczi,  E.J.,  Halamickova,  P.,  and  Schwartz,  L.M.,  Multi-Scale  Modelling 
of  the  Diffusivity  of  Mortar  and  Concrete  , RILEM  International  Conference  on  Chloride  Penetration 
into  Concrete,  1995. 

Bentz,  D.P.,  Garboczi,  E.J.,  and  Lagergren,  E.S.,  Multi-Scale  Microstructural  Modeling  of  Concrete 
Diffusivity:  Identification  of  Significant  Variables,  Cement , Concrete , and  Aggregates,  20  (1),  129-139, 
1998. 

Click  here  to  view  the  database  of  experimental  diffusivities. 


Figure  22:  Form  to  specify  inputs  for  predicting  chloride  ion  diffusivity  of  a concrete  based  on 
mixture  parameters. 
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3 Example  Applications  of  the  VCCTL 

3.1  Replacement  of  Coarse  Cement  Particles  by  Inert  Fillers 

In  low  w/c  (<  0.38)  concretes,  hydration  of  the  cement  will  be  incomplete,  as  there  is  insufficient 
pore  space  within  which  the  hydration  products  can  deposit.  Thus,  in  these  concretes,  a portion 
of  the  relatively  expensive  cement  is  potentially  being  wasted  and  serving  only  as  a reinforcing 
filler  material.  Recently,  the  VCCTL  has  been  applied  to  investigating  the  potential  of  replacing  a 
portion  of  the  coarser  cement  particles  by  inert  fillers  [7].  A few  example  results  from  this  study 
will  be  presented  here. 

As  a base  case,  the  study  considered  a concrete  with  a w/c  of  0.25.  All  microstructures  were 
based  on  the  phase  composition  and  PSD  (provided  in  the  cement  images  database)  of  Cement 
135  issued  by  the  Cement  and  Concrete  Reference  Laboratory  in  January  of  2000,  and  previously 
characterized  using  the  CEMHYD3D  hydration  model  [6].  Initial  microstructures  were  created 
with  no  cement  replacement,  and  with  the  coarsest  20.5  % and  30.8  % of  the  cement  particles 
(mass  basis)  replaced  by  an  inert  filler  material.  All  three  microstructures  were  created  using  Menu 
Selection  3 of  the  VCCTL.  The  inert  particles  were  simply  added  by  changing  the  phase  ID  for 
the  appropriate  larger  diameter  particles  from  cement  to  inert  filler  in  the  input  table  shown  in 
Figure  7.  Phases  were  then  distributed  amongst  the  cement  particles  in  each  system,  according  to 
the  measured  2-D  image  values  for  Cement  135  [6],  using  Menu  Selection  7.  The  resultant  starting 
microstructures  were  then  examined  using  Menu  Selection  4,  with  representative  2-D  slices  being 
shown  in  Figures  23  and  24  for  the  systems  with  0 % and  30.8  % cement  replacement,  respectively. 
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Slice  25  from  ceml35wc025fimg  at  5X  magnification 


Black  is  water-filled  porosity 
Red  is  C3S  (silicates)  Aqua  is  C2S 

Green  is  C3A  Yellow  is  C4AF  (aluminates ) 

Greys  are  hemi hy d rate,  gypsum,  and  anhydrite 
Orange  is  C-S-H  Blue  is  CH 

Pale  green  is  other  hydration  products  and  aluminosilicates 
Magenta  is  aggregate  White  is  empty  porosity  (and  fly  ash) 

Blue-grey  is  inert  filler 

Return  to  the  main  menu 

Figure  23:  2-D  slice  from  3-D  starting  microstructure  for  Cement  135,  w/c=0.25,  with  no  cement 
replacement  by  inert  fillers. 
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Slice  50  from  cl35wc025fil30f.img  at  5X  magnification 


Black  is  water-filled  porosity 
Red  is  C3S  (silicates)  Aqua  is  C2S 


Green  is  C3A  Yellow  is  C4AF  (aluminates) 

Greys  are  hemihydrate,  gypsum,  and  anhydrite 
Orange  is  C-S-H  Blue  is  CH 

Pale  green  is  other  hydration  products  and  aluminosilicates 
Magenta  is  aggregate  White  is  empty  porosity  (and  fly  ash) 

Blue-grey  is  inert  filler 

Return  to  the  main  menu 

Figure  24:  2-D  slice  from  3-D  starting  microstructure  for  Cement  135,  w/c=0.25,  with  30.8  % 
cement  replacement  by  inert  fillers. 
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The  systems  were  then  hydrated  for  4000  dissolution/reaction  cycles  [2],  representing  about  200 
d of  real  hydration  time  under  saturated /sealed  conditions  at  25  °C  [6],  using  Menu  Selection  13. 
Sat urated /sealed  conditions  indicate  that  the  hydration  is  initially  executed  under  saturated  (water 
is  imbibed  into  the  3-D  microstructure)  conditions,  but  the  curing  is  switched  to  sealed  curing  (no 
further  water  imbibition)  when  the  water-filled  porosity  depercolates  during  the  hydration  process. 
For  this  study,  the  percolation  characteristics  of  the  water-filled  capillary  porosity  were  evaluated 
after  every  100  cycles  of  the  hydration  process.  Representative  images  of  2-D  slices  after  the 
hydration  are  provided  in  Figures  25  and  26.  In  the  latter  image,  the  hydrated  microstructure  is 
seen  to  be  uniformly  dense,  even  in  the  immediate  vicinity  of  the  inert  filler  particles,  due  to  the 
low  overall  w/s  value. 
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Slice  25  from  ceml35wc025f.img.4000.25.000  at  5X  magnification 


Black  is  water-filled  porosity 
Red  is  C3S  (silicates)  Aqua  is  C2S 


Green  is  C3A  Yellow  is  C4AF  (aluminates ) 

Greys  are  hemi hy d rate,  gypsum,  and  anhydrite 
Orange  is  C-S-H  Blue  is  CH 

Pale  green  is  other  hydration  products  and  aluminosilicates 
Magenta  is  aggregate  White  is  empty  porosity  (and  fly  ash) 

Blue-grey  is  inert  filler 

Return  to  the  main  menu 

Figure  25:  2-D  slice  from  3-D  hydrated  microstructure  for  Cement  135,  w/c=0.25,  with  no  cement 
replacement  by  inert  fillers. 
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Slice  50  from  cl35wc025fil30f.img.4000.25.000  at  5X 

magnification 


Black  is  water-filled  porosity 
Red  is  C3S  (silicates)  Aqua  is  C2S 

Green  is  C3A  Yellow  is  C4AF  (aluminates) 

Greys  are  hemihydrate,  gypsum,  and  anhydrite 
Orange  is  C-S-H  Blue  is  CH 

Pale  green  is  other  hydration  products  and  aluminosilicates 
Magenta  is  aggregate  White  is  empty  porosity  (and  fly  ash) 

Blue-grey  is  inert  filler 


Return  to  the  main  menu 

Figure  26:  2-D  slice  from  3-D  hydrated  microstructure  for  Cement  135,  w/c=0.25,  with  30.8  % 
cement  replacement  by  inert  fillers. 
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The  predicted  degrees  of  hydration  and  compressive  strength  developments  are  compared  amongst 
the  three  systems  in  Figures  27  and  28,  which  provide  comparisons  of  the  base  system  (no  replace- 
ment) vs.  the  systems  with  20.5  % and  30.8  % replacement,  respectively.  Compressive  strengths 
for  ASTM  C109  [14]  mortar  cubes  have  been  estimated  using  Power’s  gel-space  ratio  theory  [11], 
assuming  that  the  inert  filler  particles  and  unreacted  cement  cores  contribute  equally  to  the  com- 
pressive strengths  of  the  mortar  cubes.  The  systems  with  the  inert  filler  replacements  are  observed 
to  attain  larger  degrees  of  hydration,  due  to  their  effectively  higher  w/c  ratio.  While  they  have  less 
cement  reacting  overall,  this  increase  in  hydration  rate  results  in  a predicted  compressive  strength 
development  that  is  inferior  to  but  on  the  same  order  as  the  system  with  no  cement  replacement. 
The  actual  decreases  in  predicted  compressive  strengths  are  plotted  in  Figure  29.  For  both  replace- 
ment levels,  the  difference  in  compressive  strength  is  observed  to  reach  a maximum  after  about 
14  d and  then  decay  towards  zero.  While  the  > 10  MPa  maximum  predicted  strength  loss  with 
the  30.8  % replacement  level  might  be  deemed  unacceptable,  the  6 MPa  maximum  loss  with  the 
20.5  % replacement  level  is  most  likely  a reasonable  tradeoff  for  the  cost  savings  that  would  result 
from  replacing  relatively  expensive  cement  by  a hopefully  inexpensive  basically  inert  filler  (such  as 
limestone). 


Figure  27:  Predicted  degree  of  hydration  and  compressive  strength  development  for  base  w/c=0.25 
system  and  system  with  coarsest  20.5  % of  cement  particles  replaced  by  inert  filler.  Solid  line  is 
original  system  and  dotted  line  is  system  with  20.5  % replacement. 
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Figure  28:  Predicted  degree  of  hydration  and  compressive  strength  development  for  base  w/c=0.25 
system  and  system  with  coarsest  30.8  % of  cement  particles  replaced  by  inert  filler.  Solid  line  is 
original  system  and  dotted  line  is  system  with  30.8  % replacement. 
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Reduction  in  strength  (MPa) 


Time  (h) 

Figure  29:  Predicted  reduction  in  compressive  strength  due  to  coarse  cement  particle  replacement 
for  Cement  135  with  w/s=0.25.  Solid  line  is  reduction  with  20.5  % replacement  level  and  dotted 
line  is  reduction  with  replacement  level  of  30.8  %. 
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3.2  Seeding  of  the  Hydration  with  Fine  C—S—H  particles 

The  VCCTL  can  also  be  used  to  explore  unconventional  means  of  modifying  cement  hydration 
kinetics  and  microstructure  development.  For  example,  the  use  of  microfine  C—S—H  particles  to 
seed  the  hydration  of  cement  has  recently  been  investigated  [15].  Here,  a simple  comparison  can 
be  made  between  the  predicted  hydration  and  strength  development  of  systems  with  and  without 
the  C—S—H  seeding.  Once  again,  Cement  135  is  used,  this  time  with  a w/c  of  0.40.  The  base 
system  is  compared  to  one  with  an  addition  of  0.0026  g C—S—H/ g cement,  a very  low  addition  rate. 
Specifically,  this  required  the  addition  of  only  1717  one-pixels  particles  of  C—S—H  to  the  starting 
microstructure  when  specifying  the  input  for  Menu  Selection  13,  via  the  entry  just  beneath  the 
“Inert  filler’  row  in  the  middle  of  Figure  16.  The  phase  ID  used  for  C—S—H  is  13,  as  indicated  in 
the  table  in  Figure  11. 

Both  systems  were  then  hydrated  under  saturated /sealed  conditions  at  25  °C  for  1600  cycles 
(or  about  32  d)  of  hydration  using  Menu  Selection  13.  The  predicted  degrees  of  hydration  and 
compressive  strength  developments  are  shown  in  Figure  30.  Because  the  C—S—H  concentration 
affects  the  induction  period  in  the  CEMHYD3D  model  [1],  addition  of  the  C—S—H  seeds  is  seen 
to  shorten  the  induction  period  and  increase  the  predicted  degrees  of  hydration  and  compressive 
strengths  at  all  times.  After  4 days,  the  projected  increase  in  compressive  strength  due  to  the 
C—S—H  seeding  is  on  the  order  of  5 MPa,  while  after  32  days  it  has  decreased  somewhat  to  about  3 
MPa.  Particularly,  the  hydration  of  the  cement’s  silicate  phases  are  increased  by  the  seeding  with 
microfine  C—S—H  particles.  After  32  d,  the  degree  of  hydration  of  the  C%S  is  increased  from  0.8876 
to  0.9088,  while  that  of  the  C^S  is  increased  from  0.6198  to  0.6650.  These  degrees  of  hydration 
of  specific  phases  were  easily  determined  by  using  Menu  Selection  9 to  compute  the  phase  counts 
before  and  after  hydration  for  each  of  the  two  systems.  This  simulation  study  is  very  preliminary  in 
nature,  but  does  suggest  that  the  benefits  of  increasing  silicate  hydration  by  the  addition  of  “seed” 
C—S—H  particles  could  be  significant. 
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Degree  of  hydration 


Cement  135  original  w/c  = 0.40 
Solid  line-  original  Dotted-  C-S-H  seeded 


Figure  30:  Predicted  increase  in  hydration  and  compressive  strength  due  to  “seeding”  microstruc- 
ture with  1 fim  C-S-H  particles. 
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4 Summary 


A detailed  description  of  the  web-based  menu-driven  interface  to  Version  1.0  of  the  VCCTL  has 
been  provided.  The  interface  consists  of  a series  of  fill-in  HTML  forms  that  allow  the  user  to  specify 
inputs  for,  and  view  results  of,  executing  the  various  component  programs  of  the  VCCTL.  Users  are 
notified  of  the  completion  of  program  execution  and  in  one  case  sent  a copy  of  the  analysis  results 
by  e-mail.  The  system  has  been  designed  to  bring  the  VCCTL  to  the  user’s  desktop  computer, 
where  their  Internet  browser  provides  them  with  a powerful  connection  to  the  world  of  virtual 
design,  analysis,  and  optimization  of  cement-based  materials. 

The  applications  presented  above  hopefully  provide  a brief  insight  into  the  potential  power  of  the 
VCCTL  for  conducting  research  and  development  of  cement-based  materials.  The  microstructure 
models  underlying  the  VCCTL  will  continue  to  be  enhanced  and  new  features,  such  as  the  prediction 
of  elastic  and  rheological  properties,  are  currently  under  development. 
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